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Decoherence in open Majorana systems 
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Coupling to a thermal bath leads to decoherence of stored quantum information. For a system 
of Gaussian fermions, the fermionic analog of linear or Gaussian optics, these dynamics can be 
elegantly and efficiently described by evolution of the system’s covariance matrix. Taking both 
system and bath to be Gaussian fermionic, we observe that decoherence occurs at a rate that is 
independent of the bath temperature. Furthermore, we also consider a weak coupling regime where 
the dynamics are Markovian. We present a microscopic derivation of Markovian master equations 
entirely in the language of covariance matrices, where temperature independence remains manifest. 
This is radically different from behaviour seen in other scenarios, such as when fermions interact 
with a bosonic bath. Our analysis applies to many Majorana fermion systems that have been 
heralded as very robust, topologically protected, qubits. In these systems, it has been claimed 
that thermal decoherence can be exponentially suppressed by reducing temperature, but we find 
Gaussian decoherence cannot be cooled away. 


Thermalization through interaction with an external 
bath is one of the principal mechanisms by which quan¬ 
tum systems lose information. In quantum technologies, 
rapid thermalisation destroys their advantage over clas¬ 
sical counterparts. By better understanding these pro¬ 
cesses, one hopes to identify and engineer physical sys¬ 
tems that act as more robust stores of quantum infor¬ 
mation. In topologically ordered systems, information 
is stored non-locally within the degenerate ground space 
of some large many-body system. The primary benefit 
of topology is robustness against random adiabatic fluc¬ 
tuations in the system Hamiltonian. Damage from such 
noise is exponentially suppressed with system size. Topo¬ 
logical systems also have an energy gap A between the 
degenerate ground space and excited states, and are said 
to be protected by the gap against thermal excitations. 
A common claim [1] is that thermal processes occur at 
a rate , which is sometimes called the Arrhenius 

law. The bold conclusion is that topology can exponen¬ 
tially eliminate noise merely by increasing system size 
and decreasing temperature. 

Of all topological systems, Majorana zero modes have 
attracted the most attention. It was theorized that a 
so-called Kitaev wire supports Majorana zero modes at 
edges, which could be realised in simple solid state het- 
rostructures [2], for example a nanowire coupled to a 
conventional s-wave superconductor [3] . This drove a se¬ 
ries of experiments, eventually leading to observations 
of Majorana edge modes Hi]. Beyond topological ro¬ 
bustness, Majorana zero-modes also possess the braiding 
statistics of non-Abelian Ising anyons. Though insuf¬ 
ficient for direct quantum computation, braiding Ising 
anyons can demonstrate nonlocality, teleportation and 
superdense coding [7]. Furthermore, Ising anyon braid¬ 
ing can be promoted to full quantum computing when 
supplemented with some nontopological (noisy) opera¬ 
tions [HUH. 
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The physics of these Majorana systems is especially 
tractable as their Hamiltonians are quadratic in fermion 
creation and annihilation operators. We say such a sys¬ 
tem is Gaussian, or quasifree fermionic, in analogy with 
Gaussian linear optics. Gaussian states can be described 
purely in terms of the expectation value of quadratic ob¬ 
servables, which are captured by a covariance matrix. 
Furthermore, some dissipative processes can be described 
within this powerful covariance matrix formalism nnHis], 
and allow single fermions to hop between system and 
bath via a^a^- Single fermion hopping violates conser¬ 
vation of fermion parity in the system, which is other¬ 
wise respected by unitary evolution. It is a toxic pro¬ 
cess that can cause errors without creating excitations, 
circumventing arguments that energy penalties suppress 
thermal processes to a rate . In particular, Ma¬ 

jorana modes in the Kitaev wire (see Fig. have been 
shown to decohere due to fermion hopping at rates inde¬ 
pendently of system size or the system gap naim, and 
we will review these results in detail. This article consid¬ 
ers all Gaussian fermionic systems, not just the Kitave 
wire, and how they decohere as a function of tempera¬ 
ture. A single fermion appearing in the system will have a 
partner appear in the environment, and so perhaps there 
is hope that a gapped bath Hamiltonian will provide an 
energy penalty inhibiting these processes. 

We present a very general, yet simple, argument that 
decoherence is independent of temperature, assuming 
only that the system-bath is governed by a Gaussian 
Hamiltonian. We extend this argument by providing a 
microscopic derivation of a master equation in the weak 
coupling regime, and again observe temperature inde¬ 
pendent decoherence. When fermions couple to bosonic 
baths, or through quartic fermion-fermion interactions 
one would instead find a non-trivial tem¬ 
perature dependence. Any real physical system will ex¬ 
perience noise from multiple sources, mostly with tem¬ 
perature dependent rates. However, fermionic hopping 
presents a constant background noise that cannot be sup¬ 
pressed through cooling. This adds to a growing body 
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FIG. 1. A chain of cool Majorana fermions weakly coupled to a large thermal 2-dimensional bath of Majorana fermions. This 
provides an example of the general paradigm we work within, though our results apply to all Gaussian systems. 


of work [T 9 H 22 ] that shows the outlook for Majorana 
fermions makes them less promising quantum memories 
than initially supposed. Our conclusions can be avoided 
by going beyond Gaussian fermions, for instance by mak¬ 
ing use of complex many-body interactions used in the 
passive quantum memories reviewed in Ref. [23]. We 
discuss how our results demonstrate a break down of the 
Arrhenius law, whilst still satisfying a notion of detailed 
balance. Decoherence of two-level systems, such as spins, 
has been studied when they couple to spin or fermion 
baths [111 [25| where a temperature dependence is ob¬ 
served but at low temperatures relaxation rates plateau, 
similarly breaking the Arrhenius law. Throughout, we 
use the phrase thermalize as synonymous with equili¬ 
brate or approach steady state. The reader should not 
infer that the steady state is the thermal Gibbs distribu¬ 
tion with respect the system Hamiltonian and ambient 
temperature, as may not be the case. 


1 Covariance Matrix Formalism 

Here we present standard techniques for working with 
Gaussian fermions [TOHIb] . and use them to show that 
decoherence is independent of bath temperature. Relax¬ 
ation of Gaussian fermionic open systems has be stud¬ 
ied in detail (see e.g. Refs. [Ml HU, but these did not 
include a model of the bath as also composed of Gaus¬ 
sian fermions. The first step is to map n Dirac fermions 
(e.g. electrons) with annihilation and creation operators 
into 2n Majorana operators 

C2n —1 ~ T ~ ^n)* (1) 

They are still fermionic in anti-commutation {cj^Ck} = 
but differ from Dirac fermions in that they satisfy 
cj = Cj and Cj = 1. For any quantum state p, the co- 
variance matrix has real elements composed of second 
moments 

^j,k = - CkCj)p]. (2) 

Due to fermion anticommutation, the covariance ma¬ 
trix is skew-symmetric F^ = —F. Because of conser¬ 
vation of fermion parity, first moments always vanish 


tr(cjp) = 0. For Gaussian states, expectation values of 
higher moments are determined by the covariance matrix 
via Wick’s theorem. Likewise a quadratic Hamiltonian H 
is described by a matrix H so that 

Hj^kCjCk- (3) 

iA 

Again, H must be real for the Hamiltonian to be Her- 
mitian, and furthermore H can always be chosen skew- 
symmetric H = For a closed quantum system 

evolving unitarily, the covariance matrix evolves accord¬ 
ing to 


^ = (4) 

where [•, •] is the commutator. For a time independent 
Hamiltonian, this results in 

T{t) = e^*r(0)e^"'*. (5) 


A joint system-bath covariance matrix has the form 


F = 


rc Tb T 


(6) 


where F^ and F b represent, respectively, the system and 
bath covariance matrices, and F^ records system-bath 
correlations. In other words, given a state p with co- 
variance matrix F, tracing out the bath gives a reduced 
density matrix Itb^p) with covariance matrix F^. We 
define {• • • }b to denote this process of reducing the co- 
variance matrix, so {F}^ = F^. In general, the reduced 
covariance matrix of an open quantum system will be 

Ts{t) = {e^‘r(0)e""‘}B. (7) 


For an uncorrelated system F^ = 0, the covariance ma¬ 
trix has a direct sum form F = F^ 0 F^ • The direct 
sum is linear, so that uncorrelated states have the form 
F = (F^ 0 0 ) 0 (0 0 Fb), where 0 denotes an all zero 
matrices. The covariance reduction {• • • is also linear, 
and so uncorrelated states evolve to 

^s{t) = {e^*(rB©0)e^"*}B + {e^*(0©rB)e^"*}B- (8) 
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Notice that the first term is independent of the bath vari¬ 
ables such as temperature, and can be more compactly 
written as 

{e^‘(rs © = D{t)rsD{tf, (9) 

where D{t) := We are interested in the rate 

of decoherence. How quickly will two states become in¬ 
distinguishable? Consider two different initial covariance 
matrices 0 and F^ = F|^ 0 F^, e.g. describ¬ 

ing logical encodings of qubit states. The time evolved 
difference between these covariance matrices is 

5{t) := r^(t) - r§(t) = D{t){T>’s - r^s)D{tf- (10) 

We observe that this is entirely independent of the bath 
temperature. As S{t) 0 , the states becomes indistin¬ 
guishable. Using 11 • • • 11 to denote the operator norm (the 
largest singular value) of a matrix, it is straightforward 
to show 


tr[iCjCk{p - p)] < \\S{t)\\, (11) 

where Cj and are any pair of anti-commuting Majorana 
operators. Therefore, small ||(^(t)|| entails low probabil¬ 
ity of distinguishing p and p through a single Gaussian 
measurement. We show later that this statement can be 
extended to completely general measurements. The op¬ 
erator norm is submultiplicative and transpose invariant 
so that 

||,5(i)||<||D(i)ini(rg-r§)||<2||ii(i)|^ (12) 

with smaller ||T)(t)|| entailing more decoherence. We 
have used ||(FJ^ — FJ^)|| < 2 to present an upperbound 
that is also independent of the initial state. Without 
system-bath interactions D{t) = is unitary so that 
||I)(t)|| = 1, but interactions lead to dissipation and 
||I)(t)|| < 1. Under very general conditions we have 
determined that Gaussian decoherence occurs, quite re¬ 
markably, independently of the bath temperature. In 
some instances ||T^(t)|| may decrease with time, only to 
revive later. However, for a sufficiently complex bath we 
expect Markovian behavior lead to exponentially fast de¬ 
coherence ||T^(t)|| = for some A. Next we introduce 
the formalism of Gaussian fermionic master equations, 
and then proceed to perform a microscopic derivation 
assuming weak-coupling. In such derivations various ap¬ 
proximations are made, yet we find they respect temper¬ 
ature invariance. 


2 Gaussian fermionic master equations 


Here we review Gaussian fermionic master equations 
following Refs. nni El ng. The dynamical equation for 
such a system is 


dV{t) 


XT{t) + T{t)X'^ + Y, 


(13) 


where, for dissipative dynamics, X is not necessarily 
skew-symmetric. In general, X = —H — P where H rep¬ 
resents the unitary component (and is skew symmetric) 
and P represents the dissipative part (and is symmetric 
and positive). Typically, dissipative systems will have 
one steady state Tss that satisfies XTss + F^gX^ = —F, 
and then 


m = e^*(r(0) - + r,,. (i4) 

These very general dynamical equations are the covari¬ 
ance matrix representation of a class of Lindblad master 
equations of the form 

^=i[H,p{t)]+J2{2L,pit)Ll-{LlL„p{t)}), (15) 

fl 

where the Hamiltonian H must be of the quadratic form 
introduced in Eq. and the Lindblad operators are 
linear in Majorana operators 

Lfx = (lb) 

3 

One can prove nmg that such a master equation gives 
rise to a Gaussian quantum channel with 

X =-iF-(M0M*), (17) 

F = 4i(M0M*), (18) 

where M has elements Mj^k = matrices 

X and Y will always be real so that X^ = X^ and = 

Y^. In the literature, there is some discussion of how 
such systems decohere Majorana qubits. However, prior 
work has centered on only the Kitave wire and has not 
considered how variables in Lindblad equations depend 
on the underlying system parameters. 


3 Rates of thermalisation 

We quantify the decoherence by the operator norm of 
S{t) := r^(l) -r§(i). In general, ||(5(1)|| < 2||Z)(i)||2 
and from the previous section we see that D{t) = 
Theorem IX.3.1 of Bhatia [28] shows that for any X, we 
have ||e^^|| < as previously employed in 

this setting by Bravyi-Konig m- Using the shorthand 
P := (X 0 X’f')/2, this entails 

\\Sit)\\ <2||e-^‘||2 =2e-2^^‘. (19) 

We know Ap is non-negative, so provided P does not have 
any zero eigenvalues there will be exponentially rapid de¬ 
coherence. It is well known that decoherence rates are 
governed by the spectral properties of P. Our contribu¬ 
tion is investigation of the dependence of this spectrum 
on microscopic factors. 

The above statement holds for any pair of Gaussian 
states and any Gaussian channel. We are also interested 


dt 
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in the more specific scenario where the Gaussian states lie 
in the degenerate groundspace of a physical Hamiltonian. 
Eigenstates of this Hamiltonian can be simultaneously 
diagonalised, and similarly their covariance matrices can 
be simultaneously brought into Williamson normal form. 
There exists a orthogonal matrix O, such that for every 
pure Gaussian that is a eigenstate of the Hamiltonian 
we have 


rV' = o 07/p 


( 20 ) 


two states can always be captured by an observable M 
with eigenvalues ± 1 , with an average success probability 

Pi'= ^(l + tr[M(p-p)]). (25) 

It is well known that Pr ^ 2 !!^ where the trace 

norm is ||^||tr •= tr[V A^A]. Therefore, we aim to show 
convergence in 1 -norm. We again compare two initial 
pure states encoding a qubit into 4 Majorana modes, and 
find that the time evolved density matrices states p{t) 
and p{t) satisfies 


where 7 ^ G {± 1 } distinguish different eigenstates, and 


P = 


0 1 
-1 0 


( 21 ) 


||p-p||tr < 2e 


-2Xp^t 


(26) 


This follows quickly from Eq. (24) as we show in App.[A| 


Eor some values of j, the numbers tell us whether 
there is an excitation present. However, for degenerate 
Hamiltonians there are a set of degeneracy indices G, 
such that can vary without creating excitations for 
all j G G. Therefore, the covariance matrices of ground- 
states break up into two blocks T^ = 0[r^ 0 Te\0^ 
where Te is the same for all groundstates, and 

rS = ©7*(i- ( 22 ) 

jeG 

Let us consider two encoded ground states \ 2 p) and |^). 
We deduce that 5 = = 0[{r% - © 0]O'^ 

where 0 is a zero matrix. 

The matrix representing Hamiltonian dynamics has 
the same block structure as the covariance matrices, so 
that H = 0[Hg 0 He]0^. When such a system is ex¬ 
posed to an environment, its dynamics are dictated by 
some matrix X. In many situations (including weakly 
coupled Markovian systems), X will obtain the same 
block structure as iL, so that X = 0[Hg 0 Xe]0^. It 
follows that 

5{t) := © 0]O^. 

(23) 

Defining Pq = {Xg + Xq)/2, and using the same argu- 
ments as earlier we find 

\\Sit)\\< - r^)|| < (24) 

Therefore, the decoherence of encoded ground states is 
governed by the spectrum of Pq = —{Xg + Xq)/2. 

The above arguments tell us that two initial Gaussian 
states undergoing Markovian dynamics will exponentially 
converge towards having identical covariance matrices. 
Therefore, the probability of distinguishing these states 
through a single Gaussian measurement decreases expo¬ 
nentially in time. However, a non-Gaussian measurement 
or multiple Gaussian measurements could prove more 
successful. In general, any strategy for distinguishing 


4 Derivation of master equation 


In this section we present a weak-coupling derivation 
of a Gaussian master equation in the covariance matrix 
formalism. Many of the steps directly mirror those made 
in a textbook density matrix derivation (see e.g. Ghap 3 
of Ref. [29]). We assume both the system, the heat bath 
and their interaction is entirely Gaussian. In addition, we 
make the usual assumptions involved in deducing master 
equations, notably that system-bath coupling is weak and 
that the system-bath are effectively uncorrelated at all 
times. The whole system-bath dynamics are described 
by a Hamiltonian with block matrix structure 


H = 


Hs -Hj\ 

Hi Hb 


(27) 


where Hs and Hb represent, respectively, the system and 
bath Hamiltonians and satisfy = —H^ ioi x = S^B. 
The interaction is represented by Hj a real-valued, not 
necessarily square, matrix. The initial (t = 0) system- 
bath convariance matrix has the form 

r(o) = ) . (28) 


Before proceeding we shift to an interaction picture, 
defining rint(t) := follows 

that 

= [Pint(t),/lint (i)], (29) 

where 


Hintit) 


„{Hs<SHB)t ( 0 -Hi \ (Hs^HbY 

\Hi Q 


(30) 


This simplifies to 


//int(/) 


0 -Hjit) \ 
Hi{t) 0 ) ’ 


(31) 
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where Hj(t) = Once in the interaction 

picture, we integrate over time to find 

TintW = hint (o)+ Ar int {s),H,^,{s)]ds. (32) 

Jo 

Therefore the time derivative is 

Jo 

We are only interested in the system covariance ma¬ 
trix, which is the covariance reduction of fhe 

above expression. It is straightforward to verify 
{[rint( 0 ),iJint(i)]}B = 0 , SO 

= j\ii^intis),HUs)],HUt)]}Bds. (33) 

Next, we assume the coupling is weak and that the sys¬ 
tem stays uncorrelated from the bath. Formally, this 
entails that rint(<s) ^ {rint(A}B 0 and also that 
— s) and the integral is extended to in¬ 
finity. Such assumptions directly mirror those made on 
the level of Hilbert spaces and result in the expression 


spectrum may become continuous and the integrand may 
vanish in the large s limit. Furthermore, to yield Marko¬ 
vian dynamics the resulting X and Y must be time inde¬ 
pendent. Before proceeding we can already observe that 
all F B dependence has vanished from X. 

Presently, the matrix X still carries an overt time de¬ 
pendence, which can be removed by making the secular 
approximation (SA). First we recall the explicit time de¬ 
pendence, Hi{t) = so that 

POO 

Jo 

POO 

= / (37) 

Jo 

We proceed by noting that Hs is real and skew- 
Hermitian, so it has imaginary eigenvalues iujj, eigen¬ 
vectors |j), and a diagonal form 

Hs = iJ2coj\j){j\. (38) 

3 

This entails 

/ CX3 

3,!^ - 

where 


At = [[rint(i), ffint(i - s)], ffint(i)]}Brfs, (34) 

where have made use of the shorthand T{t) := {Tint W}b. 
Next we may use our knowledge of the block structure of 
the covariance matrices to evaluate the commutator, and 
find 


{[[ri„t(t),77int(i-s)],Fi„t(t)]}B = 

- Hj- s)f{t) - {t - s)Hi{t) 

+ Hj-s) + Hj(t - s)TBHi{t). 


Combing this with Eq. 
match Eq. (13) we have 


3, and collecting terms to 


(35) 


=Xt + TX'^ + Y, 
at 


where 


X 


POO 

- s)ds, (36) 

Jo 

POO 

Y= / Hf {t)rBHi{t - s) + Hf {t - s)TBHi{t)ds. 
Jo 

We have succeeded in deriving a form of a Gaussian quan¬ 
tum channel. Though to make these equations meaning¬ 
ful we require that the integrals converge to finite values. 
Eor finite size matrices the integrands will be periodic 
functions and typically do not converge to a finite value. 
Whereas, in the limit of infinite matrices the eigenvalue 


fj,kis) = {j\Hje-^-^Hi\k). (39) 

The SA asserts that terms with rapidly oscillating phases 
^i{ujj-ujk)t neglected, except of course when ujj — 

cj/c = 0. This is valid at times longer than the reciprocal 
of the smallest energy gaps, t ^ [min^^.^^^ \ujj — • 

Eor now, we assume this to be true, but later we will 
show that a much weaker energy gap condition entails 
many of the same features. Making the SA leads to: 

POO 

^ = -E E / ^'"’^^fi,k{s)ds (40) 

7 Jo 

j k]UJk=ujj 

We see that SA has removed any dependence on t making 
time evolution Markovian. Eurthermore, the SA forces 
X to commute with Hs, and so X has the same block- 
diagonal structure as 

All matrices can be decomposed into X = —H — P 
where = —H and = P. Performing just such a 
decomposition of X we can show, via Bochner’s theorem, 
that the Hermitian part P has eigenvalues that are real 
and nonnegative (see App.[^. Eurthermore, both matri¬ 
ces have real-valued elements, so and P^ = P^, 

which entails that H has purely imaginary eigenvalues, 
just as expected. In Sec. we saw that decoherence rates 
are governed by P. Recall that Pq is the restriction of 
P to the kernel of (equivalently the groundspace of 
the associated Hamiltonian, assuming Eq = 0), and the 
decoherence rates between groundstates are governed by 
the spectrum of Pq- Eurthermore, this restricted Pq ma¬ 
trix naturally emerges when one considers a relaxed SA 
assumption. 
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Recall that the validity of SA required that all en¬ 
ergy gaps are large compared to a relevant time scale. 
Many interesting topological systems have a degenerate 
groundspace with a large gap to the first excited state, 
but then the spectrum of excitations will be dense or even 
a continuum in the large system limit. This means that 
SA cannot be used to eliminate transitions between dif¬ 
ferent excited states. However, provided the groundspace 
is gapped from excitations, we will have a limited ap¬ 
plication of SA that decouples X into 0{Xg 0 Xe)0^ 
and with Xq describing the dynamics within the ground 
space. Although SA will not apply to the dynamics Xe 
of the excitations, we saw in Sec. that decoherence in 
the groundspace is governed by only Xq- In particular, 
it is dictated by the largest eigenvalue of P^. 


5 Detailed balance and the Arrhenius law 

The bath temperature only influences what state we 
converge to, and not how quickly we get there. This 
conclusion is quite remarkable. So much so, that naively 
it seems to violate some basic tenet of physics. Two 
candidates are the Arrhenius law and detailed balance. 

The Arrhenius law is an empirical rule of thumb that 
has been successful in modeling chemical reactions. Re¬ 
cently, it has been suggested that it may also apply 
to quantum memories, though violations have been ob¬ 
served in various settings [24l [25l |30]. The Arrhenius 
law predicts that decoherence times scale as where 
A is the system gap. We have a gapped degenerate 
ground space, but quasiparticles from the environment 
can poison the system without creating an excitation. 
From this perspective perhaps we should consider A = 0, 
and our temperature independence to be consistent with 
the Arrhenius law. However, although we have focused 
on ground space decoherence our observation apply also 
to the dynamics of excitations with a discrete spectrum. 
That is, the rate at which excitation populations equili¬ 
brate is also temperature independent! Indeed, the Ar¬ 
rhenius law is not a universal law and we conclude that it 
is absolutely violated in the domain of Gaussian fermions. 

Detailed balance is a form of microscopic reversibil¬ 
ity. It states that at thermal equilibrium the population 
transfer is symmetric for each process. Consequently, 
the rate of transitions must depend on the temperature. 
At first this seems to imply that decoherence rates must 
be temperature dependent. Indeed, detailed balance has 
been used to study decoherence times of spin (qubit) 
systems with the toric code Hamiltonian m and cubic 
code Hamiltonian [32]. These results revealed an expo¬ 
nential temperature dependence, and so one may expect 
this feature to be generic. Though these are highly non- 
Gaussian systems. In the Gaussian setting, Temme et. 
al have rigorously proved very general bounds on 
thermalization rates. Their analysis appears to show an 
explicit temperature dependence, but closer inspection 
reveals variables that depend on the specific features of 


the bath. These variables can be set, whilst still respect¬ 
ing detailed balance, to precisely cancel all temperature 
dependence. We have seen that when the whole system- 
bath is Gaussian, this is indeed what happens. To fur¬ 
ther illustrate that temperature independence is consis¬ 
tent with detailed balance we present in App. a very 
simple classical Markov process where convergence rates 
decouple from temperature. Detailed balance has an es¬ 
teemed history going back to Boltzmann, who proclaimed 
it a key axiom of statistical mechanics and used it to great 
effect. However, there has been a recent surge of interest 
in non-equilibrium statistical mechanics, where detailed 
balance is violated, both in quantum [33l|34] and classical 
settings [35] . 

6 Comments on prior work 

Mazza et. al. [18] studied the effect of various noise 
models on the 1-dimensional Kitaev chain. Their main 
result is that various Hamiltonian perturbations will de¬ 
cohere the system, but this decoherence is suppressed 
by increasing the length of the Kitaev chain. Mazza et. 
al. conclude their paper by also discussing a more de¬ 
structive noise model, open systems dynamics (see pg.4 
of Ref. [E]). This model is a master equation with the 
Hamiltonian of the standard Kitaev wire (with chemical 
potential set to zero), and Lindblad operators 

Li,=val=V^{c2f^-i+ic2i,), (41) 

which allows for fermions hopping to the environment 
and where we use rj to parameterize the strength of the 
hopping. They give numerical plots for decoherence of 
this model, but it can also be understood analytically. 
We proceed by casting these master equations in the lan¬ 
guage of covariance matrices, and find 



Gomposing M 0 M* will cancel the imaginary parts, giv¬ 
ing 

77^ 

X = -H - {M ^ M*) = -H - ^fl, (43) 

Note that the eigenvalues of H are purely imaginary so 
that e^^ is unitary, but the dissipative component adds a 
constant real negative component. We have that D{t) = 
^xt _ Therefore, 

decoherence occurs at the rate ||D(t)|| = e~^ Here 
it is clear that system size and Hamiltonian gap play 
no role, provided rj is not a function of these variables. 
Mazza et. al. do not discuss how 77 itself might depend 
on the energy gap or on temperature. 

Budich et. al. m made similar observations. They 
considered two models where only the end of a Ki¬ 
taev wire couples to the environment. In both models 
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the system-environment coupling has the form i^int = 
^bathCi where ^bath is some operator acting on the bath. 
They also considered the standard Kitaev chain with zero 
chemical potential, so that the interaction commuted 
with the system Hamiltonian. They observed rapid deco¬ 
herence of Majorana edge modes, with no dependence on 
the energy gap. These toy models are excellent ways to 
illustrate a serious deficit in prior claims to effectiveness 
of topological protection of Majorana edge modes. How¬ 
ever, they tell us little about what to expect when the 
interaction and Hamiltonian do not commute. Further¬ 
more, one may also wish to consider much more exotic 
models, such as Majorana fermions in 2D systems or even 
higher dimensions. These gaps in previous work are now 
filled by the more general insights presented here. 

7 Conclusions and Acknowledgements 

We have seen that decoherence of Gaussian fermionic 
systems cannot be reduced by cooling. For Markovian 
dynamics, we provided a microscopic derivation of the 
master equation in the weak coupling limit, leading to 
exponentially fast decoherence. Therefore, to use Gaus¬ 
sian systems as a quantum memory they must be either 
highly non-Markovian or have minimal tunneling with 
any nearby Gaussian heat baths. Eliminating tunneling 
is potentially challenging when in any of the many popu¬ 
lar proposals for acquiring topological order through the 
proximity effect [ 31137 ]. In these proposals, electron hop¬ 
ping with an external s-wave superconductor is the mech¬ 
anism by which topological robustness is acquired. Both 
electron hopping and the superconducting Hamiltonian 
are Gaussian, and so this opens the door to temperature 
invariant decoherence. For such systems it is urgent that 
we acquire a better understanding of the proximity effect 
from an open systems perspective. 

After completing this work, the author forwarded the 
manuscript to Leonardo Mazza who in return kindly 
shared several unpublished yet interesting results j38F 
BO]. These tackle related problems of interactions with 
fermionic baths, including various numerical simulations 
of small fermionic (and bosonic) baths and numerous an¬ 
alytic insights. Particularly relevant is Sec. 7 of his PhD 
thesis [38], where Mazza makes several observations also 
made here, although does he not remark on the temper¬ 
ature independence of decoherence rates. 

We thank Pieter Kok and Keith Burnett for interest¬ 
ing discussions on Majorana fermions that lead to this 
research. We thank Michael Kastoryano, Tomaz Prozen, 
Jens Eisert and Leonardo Mazza for comments on the 
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A Trace norm convergence 

We assume initial states of the form p(0) = |'0)('0| = 
|'0 g)('0g| G Mq± and p(0) = |0)(0| G Mq± encoding dif¬ 


ferent qubit states. These initial states differ only on 
4-Majorana modes within the groundspace, and the evo¬ 
lution of the covariance matrix shows that this property 
holds at later times so that 

pit) = pgH) 0 Mq± , Pit) = pcit) 0 Mg±. (A1) 
Therefore, 

\\pit) - Kt)\\u = Wipcit) - Pait)) ® MG±\\tr , . 

= \\pGit)-p'ait)\\tr, ^ ’ 

where we have used \\A (g) H||tr = ||A||tr||5||tr and 
||M(^i||tr = 1. The Hilbert space of 4 Majorana modes 
supports one qubit in the even parity subspace and one 
qubit in the odd parity subspace. In other words, pcit) = 
Pg\'^)^Pg\^) and similarly poCt) = Us¬ 

ing ||A0 H||tr = ||A||tr + ll^lltr wc havc 

\\pit) - Kt)\\u = E (A3) 

x=0,l 

Eor a single qubit, we have ||p||tr = maxp^ 23 tr[pp] where 
the maximum is over all single qubit Hermitian unitary 
operators, such as the Pauli spin operators. In a Majo¬ 
rana encoding the Pauli spin operators, indeed all single 
qubit Hermitian unitary operators, are quadratic observ¬ 
ables. These expectation values never exceed the opera¬ 
tor norm of the S{t). Therefore, 

||pW-pWI|tr<||(5W||. (A4) 

Einally, we make use of Eq. p4| ) to arrive at Eq. ( |2^ . 

B Remarks on detailed balance 

Here we describe the concept of detailed balance for 
2-state systems. We show in this simple setting the con¬ 
cept is consistent with temperature invariant decoher¬ 
ence rates. Eurthermore, we show that Gaussian 2-mode 
Markov processes always obey this principle. Gonsider, 
a classical system with two possible states, with proba¬ 
bilities described by a Markov chain 

Eor simplicity we consider time to be in discrete steps, 
with a transition matrix 

so that v{t) = P^v. Gonserving flow of probability re¬ 
quires Pk^i 0 Pk ^2 = 1 for /c = 1, 2. 

We say tt is a stationary state of P, if Pit = tt, and 
where tt := (a, 1 — a). The process P satisfies detailed 
balance if 

Pl^2<^ = P2^l(l — Ol). (B3) 


One can think of tt as a thermal distribution, so a = 
Zex.p{—Ei/3) and (1 — a) = Zexp(—^ 2 /^), where Z = 
ex.-p{—EiP) + exp(—^ 2 /^) is the partition function. As 
usual, (3 is inverse temperature. In this thermal language, 
detailed balance entails that 


P2^1 

Pl^2 


= exp(A^), 


(B4) 


where A is the energy gap E^ — E\. It appears that 
the (ratio of) transition rates depend on the tempera¬ 
ture of the steady state, and so one might be tempted 
to conclude that convergence rates likewise depend on 
temperature. 

The conservation of probability and detailed balance 
give 3 independent linear constraints on P, out of the 4 
parameters of the matrix. Therefore, the space of valid 
matrices is 1 -dimensional and includes 


Pi = 


1 0 
0 1 


P = 

1 -L TT 


a a 
1 — a 1 — a 


(B5) 


Both matrices satisfy Eq. (B3). Furthermore, for all 
Markov chains v we have Pjv = v and = tt. The 
whole set of suitable matrices is contained in the span of 
these matrices, 


P^ = (1 -r])PiPvP7T 

_ / 770 ^ + (1 — 77 ) Tja \ (B 6 ) 

~V ^( 1 + ( 1 -^) y ’ 

with 0 < 77 < min[l/(a, 1/(1 — tr)] < 2 to ensure 
Pi^j G [0,1]. Within these limits, 77 is a free parameter. 
We consider a general initial probability distribution it 
always has the form 7; = tt + pp for some value p, where 
p = (1, — 1). It is easy to confirm Pip = p and Pt^p = 0. 
Therefore, Pv = (1 — r])v + p7r = tt +p(l — t])p and for t 
time steps this extends to 

v{t) = P^v = TT +p(l — pYp- (B7) 


This clearly shows exponentially rapid convergence to the 
equilibrium state at a speed governed by p. More pre¬ 
cisely, using any norm 11 ... 11 to measure distance we have 


one has a physical model of the encompassing system 
and can perform a microscopic derivation of the Markov 
process, and so derive p. This is exactly what we have 
performed for the case of Gaussian fermions, showing the 
analogous result of temperature independent p. Lastly, 
we remark that this entire discussion can be recast in 
continuous time by considering P to be the generator of 
a Markov process with transition matrix Q{t) = 

It is still interesting to ask if Gaussian Markov pro¬ 
cesses obey detailed balance. Let us just consider a pair 
of modes, with a 2 -by -2 covariance matrix 


r = 




(B9) 


The physical system is in one of two states (superpo¬ 
sitions are disallowed by fermion parity superselection), 
with probability p = (1 + A )/2 and 1 — p = (1 — A)/ 2 . A 
Markov process maps T ^ XTX^ + Y where Y is skew- 
symmetric. Under this process, we find some x^y such 
that X ^ xX P y. Therefore, p ^ px P ^{1 P y — x) and 
the probability transition matrix has the form 


P = f p^ p~ ^ py) \ 
V U^-^-y) U^px-y)) 


(BIO) 


In the steady state Pss = Pss^ + ^ (1 + ^ ~ ^) and so that 


_ I-xPy 
2(1-rc) 


_l-x-y 
~ 2(1 - x) ■ 


Therefore, 


Pss ^ l-x + y _ P 2^1 
1 - Pss 1 - a; - y Pi^2 


(Bll) 


SO that detailed balance is satisfied for all x, p. 


C Application of Bochner’s theorem 


\\v{t)-n\\=p\l-vm\. (B 8 ) 

The convergence speed is entirely independent of temper¬ 
ature, and only depends on the free parameter 77 . The 
only temperature dependence lies in 77 < min[l/o, 1/(1 — 
o)], since a depends on temperature. However, a G [0,1] 
and so the range 77 G [0,1] is valid at all temperatures. 
Therefore we can consider a family of Markov process 
with varying temperature and constant 77 G [0,1]. This 
family is consistent with detailed balance, but has a con¬ 
vergence rate independent of temperature. 

The convergence rate could vary with temperature if 
77 is a non-constant function of temperature. Although, 
there is no reason a piori to favour one function for 77 
over another. Certainly, many possible temperature de¬ 
pendencies are consistent with detailed balance. Unless, 


Here we show that the Hermitian part of X has strictly 
negative eigenvalues. The proof makes use of Bochner’s 
theorem, which relates properties of functions to their 
Fourier transform. To introduce this theorem, we first 
define the concept of functions of positive-type 
Definition 1 An absolutely integrable function (f : C ^ 
C is of positive type if for all sets of complex numbers 
{ci,C 2 ,...} the following summation is real-valued and 
positive 

^ ^ ^ra) ^ 0- (1) 

n,m 

Note that being of positive type is very different from a 
function taking positive values. Now we can state 
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Theorem 1 Bochner’s theorem: Let cj) he an absolutely 
integrahle funetion. The fourier transformed funetion f) 
is a real-valued positive funetion if and only if (j) is of 
positive type. 

Returning to the problem at hand, B = {X + 
is Hermitian by construction. From Eq. ( [^ we already 
have an expression for X, and considering X’*' we observe 
that 

POO 

j k,ujk=^j ^ 

p — oo 

= E E 1^)01/ e^'^^Ulki-s)ds, (2) 

j k,ujk=ujj ^ 

where we have made the change of variables s ^ —s. 
Switching the order of integration, 

^^ = -E E \k){j\ f e^‘^>^Ulki-s)ds. (3) 

j k,ujk=ujj ^ 

Next, we use that fj^j^{—s) = fkj{s), which can be seen 
from 

= {k\Hle-^BSHi\j), (4) 

and using = —Hb we have the result. Ap¬ 

plying this to our expression for and switching the 
dummy variables j ^ k gives 

^t = -E E 

j k,ujk=^j 


This differs from X in only the domain of the integral 
and so 

/ oo 

J K,UJk=UJj 

For each set of variables, j, /c, the integral is a Fourier 
transform of fj^k evaluated at so 

xt+x = -E E \j)m,k{s)ds. 

j k,ujk=(^j 

Notice, we have denoted Fourier transforms with a tilde. 
Next, we show that fj^k is of positive-type. For all 
{C1,C2,...} 

Tj^n(^'^fj,k{Cn - Cm) = (w^|w), (6) 

n,m 

where 

\w) = J2cme^^"'^dIi\k). (7) 

m 

Since {w\w) > 0, we can apply Bochner’s theorem and 
conclude that all fj^ki^k) are positive and real. If Hs is 
a nondegenerate matrix, then there would be no mul¬ 
tiplicity of eigenvectors with the same eigenvalue and 
—fk,k{^k) would represent the real-negative eigenvalues 
of X + Xi. However, for degenerate matrices there is a 
freedom of choice in the basis but we can always 

set this to be the eigenbasis of X+X^. Therefore, X+X^ 
has real negative eigenvalues. 


e"^^^ fj^k{s)ds. (5) 
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